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Abstract 

We summarize our recent results 1 1 , 2 1 for the induced 
exchange interaction due to thermal bosonic environment 
(bath) which also generates quantum noise. Our focus here 
is on the onset of the interaction. We demonstrate that the 
induced interaction can be used to manipulate and create 
entanglement over time scales sufficiently large for control- 
ling the two-qubit system for quantum computing applica- 
tions, though ultimately the noise effects will dominate. 

1 . Introduction 



Recently il was demonstrated I1I2I3I4I that two qubits sub- 
ject to common thermal bosonic environment (bath) can de- 
velop considerable entanglement. A similar result has also 
been obtained for qubits interacting via fermionic environ- 
ment |5j. Here we review our results on the derivation of 
the induced exchange interaction and quantum noise in a 
unified formulation 1 1 2|, focusing the presentation on the 
onset and development of the cross-qubit correlations due 
to the bath. 

We consider a ID channel model for the bath, motivated 
by recent experiments |6|, and allow bosons (e.g. phonons, 
photons) to propagate along a single direction with wave 
vector k and dispersion to k = c s k. More general results are 
available in |2|. Two qubits immersed in this environment 
are separated by distance — r\ = |d| such that the inter- 
action due to the wave function overlap is negligible. The 
qubits' interaction with the bosonic bath is introduced 1 1171 
as 

H SB = ]T 4*L (i) 

j=l,2 



where a J x is the standard Pauli matrix of qubit j 
and 



a k e 



ikr j 



1 and 2, 



(2) 



The total Hamiltonian is H = Hs + Hb + Hsb, where 
Hb = X)fc tJ '£ a ! a fc' Hs represents the Hamiltonian of the 
qubit system, and we set K = 1. The reduced density matrix 
that describes the dynamics of the qubit system is, then, 



given as the trace of the total density matrix over the bath 
modes, 



Ps(t)=Tr B (e-* Mt p s (0)p B e* Mt ) 



(3) 



where the initial density matrix is assumed factorized and 
consists of the system and bath parts. The latter is ps = 
e -H B /kT /TrB{e~ HlB / kT ). For large times, a more realis- 
tic model of the environment assumes rethermalization, and 
Markovian schemes are appropriate for the description of 
the dynamics Q. However, for short times the present for- 
mulation is adequate and provides a useful solvable model 
for the case of otherwise gapless qubits, Hs = 0, which we 
consider from now on. 

2 . Exact solution to the reduced density matrix 

With the assumptions outlined above, we utilized bosonic 
operator techniques 1 1 1 to derive an exact expression 



p s (t) = J2 p *Ps(0)Px>e* 

A, A' 



,(t) 



(4) 



Here the projection operator is defined as P\ = 
| Ai A2) (A1A2I, with \\j) the eigenvectors of crj. The real 
part of the exponent in © leads to decay of off-diagonal 
density-matrix elements resulting in decoherence, 

Re£ A v(i) = -^G fe (t,T) [(A' 1 -A 1 ) 2 + (A^-A 2 ) 2 



+ 2cos^i^l^(A' 1 -A 1 )(A^-A 2 ) 



(5) 



The imaginary part, yielding the induced interaction, is 



Im£ AA '(f) =J2C k (t) cos 



uj k Idl 



(A1A2 — A^A?,). 



We defined the standard "spectral" functions Q|8) 



Gt(! , r) = 2 M rf (^) coth (_g_ 



(6) 



(7) 



C fc (i) = 2^4- {uj k t- sin uj k t) . 



(8) 



To evaluate Q and (|6j, we consider the model in which the 
density of modes together with the coupling constants are 
approximated by the power-law function of the frequency 
with superimposed exponential cutoff Q, i.e., 



E 



9k\ 



dw Lo n exp(—uj/uj c ). 



(9) 



For n = 1 this corresponds to the well known Ohmic model 
0. 



3 . Induced cross-qubit interaction and noise 

One can show that if the real part of £\\'(t) were ab- 
sent, the exponential involving the imaginary part would 
yield coherent dynamics with the unitary evolution opera- 
tor exp[— i (H[ nt + F(t)) t]. The constant Hamiltonian H mt 
represents the induced interaction, 



Hint — 
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The time dependent term is given by 



F(t) = 2alola n I dtouj 



sin cut 
tot 



■ cos 



(11) 

F(t) commutes with H[ nt and therefore could be viewed 
as the initial time-dependent modification of the interac- 
tion during its onset: F(t) vanishes for large times as 
a n LO^/(LO c t) n , but note that F(0) = -H int . 
The interaction Hamiltonian (llOi is consistent with the 
results obtained \2\ within a perturbative Markovian ap- 
proach, for more general cases. In Figure 1, we plot the 
magnitude of the interaction Hamiltonian H.- m t, defined via 
-^int = ^int <J a; <J x> as a function of the qubit-qubit separa- 
tion for various n. At large distances the interaction de- 
creases as |d|~™, for even n, and |d|~™ -1 , for odd n. This 
means, for instance, that for spins (as qubits) with n = 1,2, 
the induced interaction decreases slower as compared to the 
dipole-dipole magnetic interaction; see estimates for semi- 
conductor impurity electron spins in 1 2 1 . 
The decoherence terms, 0, describe quantum noise that ul- 
timately destroys the coherent dynamics given by iij nt (and 
F(t)). To study the effect of these terms, we evaluate the 
concurrence |9| which measures the entanglement of the 
spin system and is monotonically related to the entangle- 
ment of formation |10|. For a mixed state of two qubits 
we first define the spin-flipped state, ps = P*s Vy^y' 
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Figure 1 : The magnitude of the induced interaction for the 
Ohmic (n = 1) and super-Ohmic (n > 1) bath models as a 
function of qubits' separation. 
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Figure 2: Concurrence as a function of time for various dis- 
tances between the qubits. The right-bottom inset demon- 
strates the topology of the concurrence in distance-time 
plane. The parameters are ot\ = 1/20, ksT/Lu c = 1/20, 
n = 1. The top inset shows dynamics for different temper- 
atures: 80k B T/uj c = 1, 2, 3, 4, 5, 6, 7, 8. 



and then the Hermitian operator R = yS/psAs^/Ps. with 
eigenvalues ^=1,2,3,4. The concurrence is then given |9| by 

C (ps) = max{0, 2 maxAi — Ai — A 2 — A3 — A4}. (12) 



In Figure 2, we plot the concurrence as a function of time 
and the qubit-qubit separation, for the (initially unentan- 
gled) state |ff), and n = 1. The bath-mediated interac- 
tion between the qubits creates entanglement, which oscil- 
lates according to the magnitude of H mt . The same bath 
also damps the oscillations destroying the entanglement for 
larger times. The decay rate of the envelope is proportional 
to the temperature, as shown in the inset of Figure 2. For 
the corresponding dynamics of the density matrix elements 
see Section 5. 



Figure 3: Initial correction to the induced interaction vs. 
time, for various distances between the qubits: o; c |d|/c s = 
0, 1, . . . , 10. The Ohmic (n = 1) case is shown. 



4 . Onset of the interaction term 

Let us now investigate in greater detail the onset of the in- 
duced interaction the time-dependence of which is given 
by F(t). In Figure 3, we plot the magnitude defined via 
F(t) = T(t)al.al, as a function of time for various qubit- 
qubit separations and n = 1. The correction is initially 
non-monotonic, but decreases for larger times as mentioned 
above. The behavior for other non-Ohmic regimes is ini- 
tially more complicated, however the large time behavior is 
similar. 

It may be instructive to consider the time dependent cor- 
rection, Hp(t), to the interaction Hamiltonian during the 
initial evolution, corresponding to F(t). Since F(t) com- 
mutes with itself at different times, it generates unitary evo- 
lution according to exp[— i / dt> 'H F (t% with H F {t) = 
d[tF(t)]/dt, 

H F (t) = al<T 2 x a n T{n) (13) 
x [u(u c \d\/c s - w c t) + u{uj c \d\/c s + w c t)], 

where u(£) = cos[n arctan(£)]/[l + £ 2 ]"/ 2 . The above ex- 
pression is a superposition of two waves propagating in op- 
posite directions. In the Ohmic case, n = 1, the shape of the 
wave is simply u(£) = 1/(1 + £ 2 ). In Figure 4, we present 
the amplitude of Hp(t), defined via Hp(t) — H.f&Io'x' as 
well as the sum of H mt and Hp(t), for n = 1. One can ob- 
serve that the "onset wave" of considerable amplitude and 
shape u(£) propagates once between the qubits, "switching 
on" the interaction. It does not affect the qubits once the 
interaction has set in. 

5 . Dynamics of the density matrix elements 

To understand the dynamics of the qubit system and its 
entanglement, let us again begin with the analysis of the 
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Figure 4: The magnitude of the time-dependent Hamilto- 
nian corresponding to the initial correction as a function of 
time and distance. The Ohmic (n = 1) case is shown. The 
inset demonstrates the onset of the cross-qubit interaction 
on the same time scale. 



coherent part in @. After the interaction, H[ nt , has set 
in, it will split the system energies into two degenerate 
pairs, E — Ei — -H[ nt and E 2 = E 3 = H{ nt . The 
wave function is then \ip(t)) = exp[—iH- mt t]\ijj(0)). For 
the initial "up-up" state, |^>(Q)) = |fT)> it develops as 
= ITT) cosftintt + |44) ismHinit. One can easily 
notice that at times t F — 7r/47ij nt , 37r/47ii nt , . . ., maxi- 
mally entangled states are obtained, while at times to = 
0, 7r/27ii nt , 7r/7ii nt , . . ., m e entanglement vanishes; these 
special times can also be seen in Figure 2. 
The bath also induces decoherence that enters via 0. The 
result for the entanglement is that the decaying envelope 
function is superimposed onto the coherent dynamics de- 
scribed above. The magnitudes of the first and subsequent 
peaks of the concurrence are determined only by this func- 
tion. As temperature increases, the envelope decays faster 
resulting in lower values of the concurrence, see the inset in 
Figure 2. 

Note also that generation of entanglement is possible only 
provided that the initial state is a superposition of the eigen- 
vectors of the induced interaction with more than one eigen- 
value (for pure initial states). For example, the initial state 
(I4T) + IT4))/ V2 in our case would only lead to the destruc- 
tion of entanglement, i.e., monotonically decreasing con- 
currence, similar to the results in 11 II . 
Since Hs — 0, there is no relaxation by energy transfer 
in the system, and the exponentials in (|4j, with l|5}, sup- 
press only the off-diagonal matrix elements, i.e., those with 
A 7^ A'. It happens, however that at large times the d- 
dependence is not important in l|3 and Re£AA'(i — > oo) 
vanishes for certain values of A ^ A'. In the basis of the 
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Figure 6: Dynamics of the off-diagonal matrix elements for 
the same system as in Figure 5. 

Figure 5: Dynamics of the occupation probabilities for the 
states ||t), |||), |jt), and ||j,). The parameters are the 
same as in Figure 1. The inset shows the structure of the 
reduced density matrix (the non-shaded entries are zeros). References 
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(14) 



in the basis of states Iff), |||), | j|), and |||). The signif- 
icance of such results, see also Q, is that in the present 
model not all the off-diagonal matrix elements are sup- 
pressed by decoherence, even though the concurrence of 
this mixed state is zero. 

The probabilities for the qubits to occupy the states |||), 
|T1)> UT)> anc l IJ--L) aie presented in Figure 5. For the initial 
state |||), only the diagonal and inverse-diagonal matrix el- 
ements are affected, and the system oscillates between the 
two states |||) and |||), as mentioned earlier in the descrip- 
tion of the coherent dynamics, while decoherence dampens 
these oscillations down. In addition, decoherence actually 
raises the other two diagonal elements to a certain level, see 
(1 1 4-1 . The dynamics of the inverse-diagonal density matrix 
elements is shown in Figure 6. 

6 . Conclusion 

To summarize, we studied the initial stages of the cross- 
qubit interaction induced by a thermal bosonic bath. It was 
shown that thermal environment can create a sufficiently 
large entanglement for quantum control, though it is erased 
for larger times. The dynamics of the entanglement and the 
density matrix elements have been investigated. 
This research was supported by the NSF under grant DMR- 
0121146. 
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